Dimer phases in quantum antiferromagnets with orbital degeneracy 
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We study and solve the ground-state problem of a microscopic model for a family of orbitally 
degenerate quantum magnets. The orbital degrees of freedom are assumed to have directional char- 
acter and are represented by static Potts-like variables. In the limit of vanishing Hund's coupling, 
the ground-state manifold of such a model is spanned by the hard-core dimer (spin singlet) cover- 
ings of the lattice. The extensive degeneracy of dimer coverings is lifted at a finite Hund's coupling 
through an order-out-of-disorder mechanism by virtual triplet excitations. The relevance of our 
results to several experimentally studied systems is discussed. 
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Isotropic quantum spin systems in dimensions higher 
than one have typicaUy magnetically long-range ordered 
ground states and thus gapless spin excitations. This 
picture may, however, fail for so-called frustrated anti- 
ferromagnets, which have extensively degenerate classi- 
cal ground states 1|. In such systems exotic quantum 
phases without long-range magnetic order can emerge as 
the true ground states. Spin systems can be frustrated 
by the geometry of the lattice or by the presence of com- 
peting interactions. In this letter, we discuss another 
scenario in which frustration of spin interactions is in- 
duced by the presence of additional electronic degrees of 
freedom, such as orbital degeneracy. A possibility of for- 
mation of orbitally driven magnetically disordered states 
has been suggested within various coupled spin-orbital 
models [lliMll. 

Here we consider a generalization of the spin-orbital 
model of Ref. Q to a wide class of lattices in arbitrary 
dimensions. We study and solve the ground-state prob- 
lem for this model. We find remarkably that orbital 
degeneracy drives the spin degrees of freedom into an 
extensively degenerate manifold of spontaneously dimer- 
ized states. We demonstrate that various types of valence 
bond crystal (VBC) states can be selected by an order- 
out-of-disorder mechanism due to virtual triplet excita- 
tions. We motivate our study by several experimentally 
discovered spin gapped systems and discuss the relevance 
of our results to these materials. 

The model- The model is defined on an arbitrary lat- 
tice with links aligned along z distinct directions. We as- 
sume that at each lattice site there are z orbitally degen- 
erate orthogonal levels such that each of them has hop- 
ing amplitudes only along the corresponding direction 
We further assume that, due to the strong on-site re- 
pulsion J7, the system is in the Mott insulating state with 
exactly one electron per lattice site. The on-site repul- 
sion also contains a small Hund's coupling Jh penalizing 
virtual processes with spin-zero state of a doubly occu- 
pied site. Such a model has experimental realizations on 
different lattices (see our discussion in the corresponding 
section). To the lowest order in t/U , where t is the hop- 



ping amplitude (assumed to be the same for all links and 
corresponding orbitals), it leads to an effective Kugel- 
Khomskii-type Hamiltonian Since second-order vir- 
tual processes locally conserve orbital Savors, the orbital 
degrees are static Potts-like variables. The effective spin- 
orbital Hamiltonian for such a system takes the form : 



H = J 



E { ■ J] - + • } (1) 



where the sum is taken over pairs of nearest-neighbor 
(NN) sites. To the leading order in Jn/U, the coupling 
constants are given by J ~ At^ /U, r) « Jh/(2C/), and 
C ^ J- The small dimensionless parameter rj character- 
izes the effective strength of the Hund's coupling. The 
full expressions for the coupling constants in terms of t, 
U , and Jh can be easily obtained from the exchange en- 
ergies given m Ref. i. For our further analysis, it will 
be important that C ^ 
coupling. 

The first term of the Hamiltonian ^ describes the an- 
tiferromagnetic (AF) coupling between NN spins and is 
active when both ends of the link are occupied by the or- 
bitals corresponding to the link direction. If denotes 
the "color" corresponding to the direction of the link ij, 



J for all values of the Hund's 



the operator Oij is defined as Oij = ^ai,aij^aj,aij where <5 
is the Kronecker symbol, and Oi — 1, z is the Potts-like 
variable for orbital flavor at the site i. The second term 
corresponds to the situation when on the bond ij only one 
orbital is of the color ay . The corresponding projector 

operator is = ^ai,ai, [1 - ^<J,,ai,\ + ^c,j,a,j [1 - ^c,,,a,,\- 

For further analysis it is convenient to rewrite the 
Hamiltonian as the sum of three terms (in the units of 
J): 



Haf + HpM + Eq , 



(2) 
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and Eq is a static energy depending on the on-site "color" 
variables, but not on the dynamic spin degrees of free- 
dom. It counts the number of matches between the colors 
of the link and one of its ends: Eq = ~( ^ ■ ■ Si^^ij ■ 

This model can be solved exactly in the case of zero 
Hund's coupling 77 = on most commonly considered 
lattices (see the following section), and a further pertur- 
bative analysis at finite ry is possible on many lattices, 
including the triangular and square lattice (see the cor- 
responding section). 

Zero Hund 's coupling - We start our analysis from the 
limit of zero Hund's coupling 77 = 0. In this limit, we 
show that the ground state is extensively degenerate and 
is described by dimer coverings of the lattice. 

The term Eq in the Hamiltonian only depends on the 
static orbital variables (colors). The antiferromagnetic 
term i^AF is active only on the links whose colors match 
those at both their ends. The active AF bonds thus form 
non-intersecting linear chains (Fig. [1] left). The longer 
is the chain, the more energy can be gained from the AF 
spin interaction. However, for each AF link we pay a pos- 
itive energy (2^ — \) > 1/4 (the second term in Haf)- 
We can prove that under these conditions the minimal 
possible AF energy is achieved when all AF chains are 
dimers. A proof follows from the variational estimate 
on the ground-state energy of the M-site AF Heisenberg 
chain with open ends: Em > j ^ with the equal- 
ity attained only at M = 2. [This estimate is, in turn, 
obtained by dividing the chain into shorter overlapping 
sub-chains of lengths two and three with exactly known 
energies E2 — —3/4 and Ez = -1.] 

Therefore, at zero Hund's coupling ry = 0, the ground- 
state coloring of vertices corresponds to a dimer covering 
of the lattice (Fig. [2 right), provided that this dimer cov- 
ering simultaneously minimizes Eq . The minimization of 
Eq is, for example, always possible for lattices where each 
site has equal valencies for all colors of outgoing links. 
This condition is satisfied for most well-known lattice 
types: triangular, square (cubic), kagome, pyrochlore, 
checkerboard, hexagonal, etc. In addition, dimer cover- 
ings minimizing Eq may also exist for other exotic lattices 
not satisfying the equal-color- valency constraint. 

The dimer coverings corresponding to the ground state 
must obey an additional "no-chain" constraint: no two 
neighboring dimers can lie on the same line (otherwise 
they form a longer AF chain which is energetically unfa- 
vorable). Even with this constraint, the entropy of such 
dimer coverings remains extensive for most commonly 
considered lattices. Thus we conclude that, on such lat- 
tices, the model is in the dimer-liquid state at zero Hund's 
coupling. 

Order out of disorder by triplet fluctuations - We now 
introduce a finite value of Hund's coupling 77 > and 
show that it lifts the extensive degeneracy of dimer cov- 
erings by triplet fluctuations. The Hund's coupling pro- 
duces a ferromagnetic term ffpM active only on links 




FIG. 1: (Color online) Left: A possible orbital pattern and 
the corresponding decoupled AF chains on the triangular lat- 
tice, in the limit of zero Hund's coupling. The numbers label 
the three bond directions and the occupied orbitals on sites. 
Right: An example of the dimer covering minimizing the en- 
ergy at zero Hund's coupling. Thick (thin) solid lines denote 
AF (FM) intra-(inter-)dimer bonds, respectively. 




FIG. 2: Examples of topologically nonequivalent linked clus- 
ters. 



with exactly one end matching the link color. Given a 
dimer covering from the previous consideration, the FM 
bonds correspond to the "legs" of dimers: adjacent links 
along the dimer direction (Fig. [1] right) . We assume that 
the Hund's coupling is small rj 1 and treat it pertur- 
batively. A convenient method to build a perturbative 
expansion is to rewrite the FM coupling between two 
ends of dimers i and j in the basis of triplet excitations 
on those dimers At each dimer, we introduce three 
triplet excitations tj ^ characterized by the spin polariza- 
tion a = x,y, z. Those triplet excitations may be treated 
as hard-core bosons, and the AF part of the Hamilto- 
nian is written (in the units of J) as Haf = a ^1 a^i-a- 
The FM part of the Hamiltonian HpM is treated as a 
perturbation and is written as the sum of the three con- 
tributions 
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ieapjmjt^'^ jljjti^^ + [i ^ j] + h.c. 
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where = ± is a phase factor depending on the orienta- 
tion of the singlet wave function on the dimer ^] . We also 
introduce A — 77/8, to simplify combinatoric coefficients. 

A convenient way to analyze perturbative corrections 
is provided by the linked-cluster expansion The 
weak inter-dimer coupling Hfm is treated as a pertur- 
bation, and the correction to the ground-state energy 
may be written as a sum over all distinct linked clus- 



ters AE = "Y^^-, NqEc- Here Nq is the number of ways 
the cluster C can be embedded in the lattice with a given 
dimer covering and Ec is the perturbative contribution 
to the energy from the cluster C (given by a power se- 
ries in A) . Examples of topologically distinct clusters are 
presented in Fig. O Note that the leading contribution 
to Ec from tree-like clusters with n FM links is of or- 
der A^" [the clusters (1), (2a), and (26) in Fig. [2], while 
loop clusters with n links [the clusters (3a), (35), (4a) 
and (45)] contribute to the order A". The actual calcu- 
lation of Ec involves lengthy combinatorial analysis of 
possible virtual processes and will be reported elsewhere. 
Here we just present the results of those calculations in 
application to the triangular and square lattices. 

The leading perturbative contribution comes from the 
cluster (1) and is given by = — 6A^(1 — 2A-f A^-|-. . .). 
However, the total number of such clusters is independent 
of the dimer covering (it is equal to the number of inter- 
dimer bonds, i.e., twice the number of dimers), therefore 
the degeneracy is not lifted to this order. 

On the triangular lattice, the degeneracy is lifted to 
the leading order by the contributions from the cluster 
types (3a) and (36) in Fig. The cluster (3a) generates 
a non-frustrated FM coupling while the cluster (35) is 
frustrated. The contributions of these two clusters differ 
in sign, -E(3a) = —-£'(36) = — 36A^, to the leading order 
in A. Thus, on the triangular lattice, the ground state 
should maximize the number of (3a) type clusters and 
minimize the number of (36) type clusters. We prove 
that such an optimization leads to a dimer crystal shown 
in Fig. [3] (left panel). It provides the maximal density 
N/ 5 (TV is the number of sites) of (3a) type clusters and 
no (36) clusters. The unit cell of such a crystal contains 
20 sites, and its degeneracy is 60-fold. A proof of crys- 
tallization may be performed by dividing the triangular 
lattice into edge-sharing hexagons composed of six tri- 
angles each. On a hexagon, we can place at most one 
triangular loop of a cluster (3a). A hexagon which has 
one of its edges shared by a triangular loop placed on an 
adjacent hexagon can not itself host a triangular loop. 
An empty hexagon can have at most two such edges (fur- 
thermore, they must be non-parallel and disconnected). 
To achieve the maximal density of triangular loops, we 
start with a three-hexagon cluster composed of an empty 
hexagon B sharing two edges with the triangular loops 
on the adjacent hexagons A and C (Fig. [3l the letters 
mark the centers of hexagons). Then the hexagon D 
must also be empty and, moreover, it may have at most 
one edge shared by a triangular loop (on the hexagon E) . 
One can show that the maximal density of (3a) clusters is 
achieved by a close-packed covering of the triangular lat- 
tice by five-hexagon clusters of the type A-E. Now start- 
ing from one such cluster and choosing different possible 
divisions of the lattice into hexagons one can verify that 
the corresponding dimer covering is unique and periodic, 
as shown in Fig. [3] (left panel). 




FIG. 3: Left: The ground state dimer covering of a triangular 
lattice. The triangular loops are explicitly shown as dark 
triangles. A possible choice of an unit cell is shown by open 
circles. Right: The ground state dimer covering of a square 
lattice. The length-four loops are shown as dark squares and 
white rectangles. 



On the square lattice, a similar analysis needs to be 
done to the order A'^, since length-three loop clusters 
are not possible. At this order, the contributions come 
from the length-two clusters (2a) and (26) and from the 
length- four loops (such as (4a) and (46)), see Fig. [2] 
(other length-four loop clusters are also possible, but not 
shown). The clusters (2a) and (26) have different contri- 
butions: E(2a) = 12A^ and £^(26) = -20A''. All length- 
four loop clusters involve unfrustrated FM couplings and 
contribute the same energy i5(4) = —120 A''. With those 
cluster energies, we can prove that the resulting dimer 
covering forms a crystal shown in Fig. [3] (right panel). 
The proof is based on the observation that the energy 
gain from length-four loop clusters is much larger than 
the scale of energy optimization between (2a) and (26) 
type clusters. Therefore, an optimal configuration should 
involve as many length-four loops, as possible (each FM 
link should belong to two such loops). Under this con- 
dition, an optimization of the energy contribution from 
clusters (2a) and (26) leads to the crystal structure shown 
in the figure. Its unit cell contains 8 lattice sites, and the 
degeneracy is 8-fold. 

On other lattices, we expect that virtual triplet pro- 
cesses also lead to a dimer crystallization at some order 
of the perturbation theory. A similar analysis may be 
performed individually for a given lattice. 

Experimental systems - The considered model and the 
derived results may be relevant for several experimentally 
studied materials. 

In the layered compound NaTi02 3, 11 1 the spin one- 
half Ti'^+ ions form a triangular lattice and have one d- 
electron in the three- fold degenerate t2g manifold {dxyj 



dxz, and dyz orbitals). The three non-equivalent bonds 
of the triangular lattice are aligned along the xy, xz, 
and yz directions, and each orbital state has a dominant 
hopping amplitude along the corresponding bond. This 
system does not show any signature of long-range mag- 
netic order at low temperatures lip|. A drop of magnetic 
susceptibility has been observed [ll|, but no dimerized 
superstructure has yet been detected. 

The compound Sr2V04 is a possible realization of the 
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model (dl) on the square lattice (formed by spin one-half 
Vanadium ions) . Each V*^ ion is surrounded by the dis- 
torted oxygen octahedra elongated in the z-direction 
and single d electron occupies a two-fold degenerate {d^z 
and dyz) level. On the square lattice in the xy plane, the 
dxz and dyz electrons hop only in the x and y directions, 
respectively. The experimental study of Sr2V04_5 
suggested an absence of magnetically ordered states down 
to zero temperature and a possible formation of singlet 
pairs in the stoichiometric compound 5 — Q. The avail- 
able experimental data are, however, inconclusive about 
the presence of a spin gap in this compound. 

The model ^ on the pyrochlore lattice describes the 
MgTi204 spinel compound |^]. The latter has a spin gap 



and a dimerized bond pattern at low temperatures [14 1. 
It has been shown in Ref. Q, within the mean- field like 
approach, that the dimer states for model ([T]) on a py- 
rochlore lattice are favorable over some other magnetic 
states. Here we provide a proof that dimer states are 
the ground states. Our mechanism of the selection of 
VBC pattern would, however, predict a different dimer 
pattern than observed experimentally in MgTi204. This 
suggests that in this compound the degeneracy is lifted 
by a magnetoelastic mechanism, as discussed in Ref. [1], 
and not by quantum fluctuations. An alternative mech- 
anism for observed dimerized pattern within a different 
(itinerant-electron) model has recently been suggested in 
Ref. [H. 

Note that our model has two different energy scales. 
At the higher energy scale (of order J), the system orga- 
nizes into classically fluctuating singlets (dimers). This 
possibly explains the spin gap in the compounds men- 
tioned above. The crystallization of dimers by triplet 
fluctuations occurs at a much lower energy scale (deter- 
mined by the rjJ term in Eq. [1]). In real systems, this 
is one of several possible mechanisms of crystallization 
(including magnetoelastic effects, residual overlaps be- 
tween orbitals in non-dominant directions, etc.) For any 
particular compound, a theoretical prediction of the crys- 
tallization pattern would require a thorough analysis of 
different ordering mechanisms. On the experimental side, 
an identification of the dimerization pattern may be dif- 
ficult in the case of large unit cells like those predicted 
in this work. 

Summary.- To summarize, we have studied a spin- 
orbital model for a family of quantum antiferromagnets 
with orbital degeneracy. In such systems, the orbital de- 
grees of freedom induce a spontaneous dimerization of 
the spins. The resulting extensively degenerate ground 
state is equivalent to a problem of constrained classical 
hard-core dimers. At a lower energy scale, the exten- 
sive degeneracy is lifted through an order-out-of-disorder 
mechanism by virtual triplet excitations and a valence- 
bond crystal is formed. At intermediate temperatures, 
above the dimer crystallization, the spin system is in a 



liquid state. Our study proposes a mechanism of the spin- 
gap formation experimentally observed in several two- 
and three-dimensional systems. An interesting possible 
extension of our model may include additional terms aris- 
ing from neglected overlaps of orbitals and leading to a 
quantum dynamics of dimers. Such an extension might 
be a better candidate than pure spi n sy stems for a re- 
alization of quantum dimer models [16l |l7| and of the 
long-sought quantum spin-liquid phase. 
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